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Abstract 

In the last years, a remarkable theoretical effort has been made in order to understand 
stability and complexity in ecological communities. The non-random structures of real 
ecological interaction networks has been recognized as one key ingredient contributing 
to the coexistence between high complexity and stability in real ecosystems. However 
most of the theoretical studies have considered communities with only one interaction 
type (either antagonistic, competitive, or mutualistic). Recently it has been proposed a 
theoretical analysis on multiple interaction types in ecological systems, concluding that: 
a) Mixture of antagonistic and mutualistic interactions stabilize the system with respect 
to the less realistic case of only one interaction type; b) Complexity, measured in terms 
of the number of interactions and the number of species, increases stability in systems 
with different types of interactions. By introducing new theoretical investigations and 
analyzing 21 empirical data sets representing mutualistic plant-pollinator networks, we 
show that that conclusions are incorrect. We will prove that the positive complexity- 
stability effect observed in systems with different kind of interactions is a mere consequence 
of a rescaling of the interaction strengths, and thus unrelated to the mixing of interaction 
types. 

Introduction 

The relationship between diversity and stability in ecosystems is one of the most debated 
issue by ecologists [j}j7]- Empirical evidences [3j[5] suggest a positive diversity - stability 
relationship, while theoretical studies challenge this point of view. The discussion started 
since publication of the pioneer work by Robert May |1| which provides a quantitative 
relation of th community stability of randomly assembled ecosystems with its number of 
species (S), connectance (C) - i.e. the ratio between actual and potential interactions 
in the ecosystem - and the characteristic interaction strength. Specifically, for random 
interactions drawn from a Gaussian distribution of zero mean and standard deviation a, 
a randomly assembled system is stable (i.e. all real part of the eigenvalues of the stability 
matrix are negative) if 

aVCS<l. (1) 
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This result rises the paradox that more complex an ecosystem is (i.e. large C ■ S), less 
stable it is [2]. This important result suggests that real networks must have some non- 
random, stabilizing structures that allow ecosystems to persist despite their complexity 
[8j[9]. The structure and the interactions types (e.g. mutualistic, antagonistic, etc..) have 
indeed a remarkable impact on the stability profiles of ecosystem community dynamics 
(7 10 ,11 . In particular it has been found that predator-prey systems are more stable with 
respect to random assembled community, while mutualistic communities are not |7||12). 
However, these results do not solve the paradox: still for structured community ecosystem 



stability decreases for increasing C or S (and fixed interaction strength) 



A possible solution to the paradox has been recently proposed 13l 14 





A theoretical 



analysis on multiple interaction types in ecological systems 13 concluded that a mixture 
of antagonistic and mutualistic interactions is more likely to stabilize the system with 
respect to the less realistic case of interactions of a single type. Moreover it was found that 
complexity, measured in terms of the connectance and the number of species, increases 
stability in systems with different types of interactions, in contrast with May's result [2]. 

In this work we show that the mixing of interaction types does not solve the paradox 
in the stability-diversity relationship. On the contrary, also for an ecosystem with hybrid 
interactions, stability decreases as complexity increases. We will first briefly summarize 
the theoretical framework and mathematical details of the community dynamics model 
used by Mougi and Kondoh 13 and variants of it. We will then show that: a) Dataset 
of real mutualistic interaction networks do no support a key assumption made in their 
framework; b) Even from a theoretical point of view a mixture of interaction types does not 
have a stabilizing effect; c) The positive complexity-stability effect recently proposed [13] 
is only a consequence of a rescaling of interaction strengths, which is unrelated to the 
mixing of interaction types. 



Materials and Methods 
Community Dynamics Model 
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The model considered by Mougi and Kondoh 13 is formulated through Lotcka-Volterra 

c 

drii 



nJri + ^a ij n j ) , (2) 



dt 

where is the abundance of species i, S is the number of species, r, is the intrinsic 
rate of growth and is the interaction coefficient between species i and species j while 
an = Zi represents the self interaction (density dependent regulation) and it is uniformly 
distributed between and 1 (z, ~ ^[0,1] )• The matrix contains all the information 
about the interactions between species. If species i has a mutualistic interaction with 



3 



species j (i ^ j), the corresponding terms in the interaction matrix will have the form 

Z^fcem(i) A ik . . 

4- ' ' 

a ji — JM e ji 



^k£m(j) Ajk 

while if species i has an antagonistic interaction with j (consider for instance i predator 
and j prey) 

a ij = f A9ij~- 



< 4 > 



_ _ , -±i 3 _ ^13 



where m(i) is the set of the species having a mutualistic interaction with i, whereas p(i) 
the species which are preys of i. Ay measures the potential preference of an interaction 
between i and j, Ja and fu are the relative strengths of the antagonistic and mutualistic 
interaction, respectively, while and quantifies the symmetry of the interactions. 
The model could also be generalized to a non-linear functional response (Holling type 
II) [II]. Not all the species interact: only a fraction C of matrix elements Ay is different 
from zero, of which a fraction pm is mutualistic, while a fraction 1 — Pm is antagonistic. 

For a given choice of i) the structure of the matrix Ay (random, cascade for antagonis- 
tic part and bipartite for mutualistic one) [7], ii) the parameters /a and fu and Hi) ran- 
domly drawn values of ey and <?y uniformly distributed between and 1, Eqs. ^ and Q 
lead to the interaction matrix a. By introducing a stationary point n* = {n\, . . . ,ns), 
whose components are randomly drawn from a uniform distribution between and 1 , and 
linearizing Eq. ^ around n* one finally obtains the stability matrix M. If all eigenvalues 
of M have negative real parts, the system is stable. Otherwise it is unstable. Mougi and 



Kondoh 13 have shown that starting from a mutualistic (antagonistic) interaction matrix 
and adding antagonistic (mutualistic) links, i.e. decreasing (increasing) pm, the stability 
increases in a non linear fashion. Moreover they also showed that for hybrid communities 
the stability increases as S and/or C increases. 

The crucial assumption made in [13], is to impose a "constant interacting effort" for 
each species, i.e. if a species is generalist and thus can positively interact with several 
different species, then the average interaction strength must be smaller than the one of 
specialist species that has only few resources. This fact should be translated mathemati- 
cally by rescaling the interactions as 

A ■ ■ A ■ ■ 

a- (5) 
Sfcem(i)u P (i) Aik SC(p M + )E(A) 
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where m(i)Up(i) is the set of the resources species of the species i, i.e. both the mutualistic 
partners or preys of i. However, in the work Mougi and Kondoh the "constant interacting 
effort" is imposed in a more stringent way: the interacting effort spent separately in 
mutualistic and in antagonistic interactions are fixed to a constant. They indeed define 
the rescaled interaction between the species % and the species j as shown in Eqs ^ 
and Q. As a consequence of this stronger assumption, the average effort spent by a 
species in mutualistic or antagonistic interactions does not depend on p M . The total 
effort spent by species i to interact with its mutualistic partners is indeed X^/emfi) a «i ~ 
/jfE(ei J -Afj )/E,Aij = }'m/2, where E(-) represents the expectation value, and we have 
used the law of large numbers and the fact and A%j are independent random variables. 
In the same way the effort spent by the species i as a predator is on average equal 
to /a/2. Therefore the authors of ref. 13 do not only assumes that, (i) "interaction 
strengths decrease with increasing resource species, due to an allocation of interacting 
effort", but also that, (ii) the total interaction strengths spent in mutualistic (antagonistic) 
interactions does not depend on the number of mutualistic (antagonistic) interactions, (i) 
and (ii) are encapsulated by the denominators in Eqs. M) and (j^J). In other words they 
are assuming that the mutualistic interacting effort spent in a ecosystem where 1% or 
99% of the links are mutualistic, is the same. 



Data Analysis 

Although at least assumption (i) might seems plausible and biologically justifiable, we 
show that this, and thus assumption (ii), is not supported by observations in empirical 
mutualistic ecosystems (see Figure fij). 

We have analyzed a dataset consisting of 21 empirical mutualistic networks [16) . All 
information on each dataset b is incorporated in the weighted adjacency matrix W b , 
whose elements gives the interaction strength between species i and j. Therefore 
the normalized total interaction strength of species i in database b is defined as s\ = 
Y.kem(i) w ik/ max kl (w b kl ), while the number of interacting partners is d\ = E fcem (i) 
(Q(x) = 1 if x = and otherwise). From these information we can analyze the empirical 
relations between s and d, i.e. the relationships between the positive number of resources 
of a species and its interaction strength. In order to minimize the effect of fluctuations due 
to intrinsic stochastic variability in ecosystem dynamics and errors in interaction strength 
detection, we have averaged this relation across all databases, i.e. 

m = E> Eg? _ (6) 

E t E£? <*«-<*) 

where S(b) is the number of species in the dataset b. 

A constant interacting effort hypothesis would suggest that s and d are anti-correlated: 
s oc 1/d. On the contrary we find that the average interaction strength and the number of 
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interactions are strongly positive correlated, giving a falsification of assumptions (z) and 
(ii). We stress that, as also Mougi and Kondoh noted [13], the relaxation of assumption 
(ii) is crucial for the validity of all their conclusions. The stabilizing effects due to the 
mixing of mutualistic and antagonistic interactions types vanish if the constant interaction 
effort hypothesis does not hold separately for the two kinds of interactions. 
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Figure 1. The interaction effort of each species is correlated to the number of 
interactions. The figure shows relation between the average interaction strength (s) and 
the number of its resource partners (cfj). We average over the normalized strength of the 
species sharing the same value of k (black points) among the databases. The error bars 
are calculated as the standard deviation of the average (cr s ). Points without error bar 
are single data points and thus they are less reliable. 



Discussion and Results 

We have shown that data do not support the assumption that interaction strengths de- 
creases with increasing resources species. We now show that even from a pure theoretical 
ground it is neither true that multiple interaction types stabilize ecosystems nor that 
increasing complexity increases stability. 
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Stability Criteria for Hybrid Communities 

We first analyze the impact of mixing mutualistic and predator-prey interactions type on 
the ecosystem stability by extending to hybrid communities the analysis recently proposed 
by Allesina and Tang [7]. In this subsection we thus go beyond the particular model 
dynamics presented above, and we use a more general approach based on . 

Consider the S x S stability matrix M for a given community dynamics hi = f(h) 
with i — 1, ...S. Mij describes the effect of the interaction between species i and j in the 
proximity of a feasible stationary point n* of the underling dynamics /: = 

^ In* 

M for hybrid predator-prey and mutualistic interactions networks can be built in the 
following way [7] . We first pick at random a pair i — j of species and we draw a random 
value p from a uniform distribution between zero and one (£7[o,i])- If p < C, the these two 
species interact, otherwise they do not interact (M^ = Mji = 0). If i — j are interacting 
species, then with probability 1-pu species i preys species j, otherwise they are mutualistic 
partners. In the former case we set ~ \Af(0,a 2 )\ and ~ — \Af(0, u 2 )\, otherwise 
~ |A/"(0, <r 2 )| and ~ |A/*(0, a 2 )\. Af(fi,a 2 ) is the normal distribution with mean 
fj, and variance a 2 , and determines the intensity of the interactions among species (with 
the notation \Af(0,a 2 )\ we mean that a random number is taken from jV(0, a 2 ) and its 
modulus | • | is taken). 

Following this simple algorithm one can build cascade predator-prey community ma- 
trices with a desired fraction pu of mutualistic interactions and then analyze the corre- 
sponding eigenvalues to study the stability for several levels of diversity of interaction type 
in the ecological community. Figure [2] shows that also for hybrid community the stability 
of the linearized dynamics described by M decreases as S and/or C increases, indepen- 
dently of pm- Moreover, as expected from previous results |7|, mutualistic interactions 
destabilize the community. 

Applying this approach to the community dynamics given by Eqs. § @ and Q it 
turns out that = —n*aij. In 13 assume that the stationary abundances n* are 



uniformly distributed between and 1 and thus stability of an hybrid community, not 
only depend on the fraction of mutualistic links, but also on their strengths. In fact, in 
the presented community dynamics framework, absolute strengths of the mutualistic and 
antagonistic interactions is defined by two parameters }m and /a- As shown in [7], in the 
large S limit there is always a real eigenvalue which is only determined by the average 
entry of the matrix 

A a = -E(M«) + (S - l)Cn ~ E(n*) [-E(z) + E(e)f M - (1 - E(g))f A j , (7) 

where \i is the average value of the out-diagonal non-zero entries of M. The second 
equality is obtained in the large S limit. Note that the right term does not depend on S 
and C, nor on p M , because of the rescaling of the interactions. By considering z iy n*, e^, 



Qij and Aij uniformly drawn between and 1, as done in 13 , this eigenvalue is greater 
than zero if /m > /a + 1- 
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Figure 2. Stability profile of the stability matrix M in the model independent 
framework presented in (7). The linearized matrix M are obtained as done in [7]. The 
profile of eigenvalues distribution in the complex plane, shows that increasing complexity 
decreases stability, i.e. the maximum real part of the eigenvalues of M increases for 
increasing S and/or C. Moreover the stability is not increased by the mixing of 
interaction types, but as expected [7] adding mutualistic links decreases stability. 
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Figure [3] shows that the stability does depend on the choice of the relative strength 
of mutualistic and antagonistic interactions, even in the case of structured interactions, 
showing that the prediction of equation ([7]) is valid in that case too. Indeed the system 
is always unstable when J'm > + These results confirm the important role that 
mutualistic interactions play in ecological networks 10 17 and poses the crucial problem 
on the estimation of these parameters from real data, a challenging but difficult task 
18lfl9 . 
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Figure 3. The relative strengths of mutualistic and antagonistic interactions strongly 
affects the stability. The figures show the eigenvalues in the complex plane for randomly 
generated matrices with S = 200 and C = 0.8 and two values of pm (while Zi, n*, Aij, 
gij and are drawn from an uniform distribution between and 1). As shown in [7], 
there is always an eigenvalue which depends only on the average interaction strength, in 
this case the eigenvalue is equal to (/m — Ia ~ l)/4 (black arrows). If /a/ > /a + 1 the 
system is, for a sufficiently large S, always unstable independently of the fraction of 
mutualistic interactions pm and the number of species S. In this case fu = 3 and 
/a — 1. If ju < /a + 1 the destabilizing effect of pm is also observed (see Figure 4). 



Scaling of the Interaction Strengths 

Figure [4] shows that if a constant interacting effort is considered (assumption (i)) the 
increased stability due to mixed interactions does not hold. This proves that the results 
obtained by Mougi and Kondoh are strongly affected by their peculiar assumption (ii). In 
other words, assuming that the total mutualistic interaction strength does not depend on 
the relative abundance of mutualistic links (assumption (ii)) is the cause of the singular 
effect on the community stability observed in (l3| . 

Finally, a positive relation between stability and complexity is also observed in |13] . 
Again, this is a consequence of the rescaling in the interaction strength and not of the 
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Figure 4. Relationship between complexity, stability and fraction of mutualistic links 
(j)m) for cascade networks [7] with Holling type I linear response. The interaction 
strengths are rescaled by summing over both mutualistic partners and preys (see Eq. |5j 
assumption (i)). Panel A shows the community stability (measured as the probability 
that the linearized matrix is stable) vs. the fraction of mutualistic links pu- Colors 
indicate different values of connectance. S is fixed to 50. Panel B shows the 
complexity-stability relationship with varying pu- Both the panels shows that mixing of 
interaction type in the community dynamics model given by Eq. ^ and based on 
assumption (i) does not increase stability. The parameters Zi, n*, Ay, gij and e y - are 
drawn from an uniform distribution between and 1, while Jm — — !• 



mixing of interaction types. In particular given assumption (ii) it can be shown (Eq. 10 
of SI in [l3]) that a system with mixed interaction type is stable if a condition of the kind 

^d<i, (8) 

is valid, where h(pM) is a suitable continuous function. This condition is exactly the 
opposite of the well known May's criterion for stability [l] given by Eq. 0. It is the 
key responsible of the Mougi and Kondoh conclusions, i.e. for hybrid ecosystem increasing 
C and/or S stabilizes the system. It is very important to note that the relation given 
by Eq. Q does not arise as a consequence of the mixture of interaction types, but it 
is instead derived from the rescaling of interactions effort at increasing resource species. 
Indeed, if the interactions are rescaled following assumption (ii), then the strength of 
of the interaction matrix elements is a decreasing function of the number of species S 
and connectance C. Figure [5] shows that also for non-hybrid ecological networks this 
stabilizing effect for increasing complexity is observed. It is exactly the imposed inverse 
correlation between SC and the interaction strengths which produces the stabilization 
effect for increasing ecosystem complexity, not a mixture of interaction types. 
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Figure 5. The rescaling of interaction strengths, and not a mixture of interaction 
types, produce a positive complexity-stability relationship. The stability is measured as 
the probability to have a positive eigenvalues. The blue curves are obtained with a 
random networks, the red curves with cascade predator-prey interactions, while the 
green ones with bipartite mutualistic matrices. Panel A and B show the dependence of 
the stability on the number of species S and the connectance C, when the interactions 
are not rescaled. As expected the stability decreases at increased complexity. Panel C 
and D are obtained with rescaled interaction strengths as in Eqs. § and Q. The 
rescaling strongly affects the stability-complexity relation. The result obtained in 13 is 
therefore not a consequence of mixing of interaction types (the green and red curves are 
obtained with purely mutualistic and cascade predator-prey interactions respectively). 
Instead it is due to the rescaling of interactions, indeed a positive complexity-stability 
relation is observable once the interactions are rescaled. 



Conclusions 

In this work we have shown that for an hybrid ecological community where a mixing 
of mutualistic and predator-prey interaction types are present, the complexity stability 
paradox is still an issue. In particular, mutualistic interactions destabilize the system as 
observed by analyzing stability profiles of structured matrices [7]. The positive relation 
between stability and complexity observed in 13 is a consequence of a particular rescal- 
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ing in the interaction strengths: by imposing that a scales inversely as the number of 
resource partners, then the stability threshold on the linearized dynamics shifts to greater 
values for greater complexity. However analysis of real mutualistic ecological networks 
support the opposite behavior: the average strength increases for increasing mutualistic 
partners. These results call for different principles beyond network structures and mixing 
of interaction types in order to understand the complexity and stability relationship in 
real ecological systems. 
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